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This report consists of two parts and appendices. The two parts are drafts of papers 
that are being prepared for publication with co-author Roger Yelle. These papers are on 
the ionosphere of Titan and heating efficiencies in the Titan thermosphere. Appended to 
the report are copies of chapters published or in press that relate to the subject of this 
project, and the grant is acknowledged, where possible, in these manuscripts. They are 
“Ion chemistry in atmospheric and astrophysical plasmas” by A. Dalgarno and J. L. Fox, 
“Dissociative recombination in planetary ionospheres” by J. L. Fox, “Hydrocarbon ions 
in the ionospheres of Titan and Jupiter” by J. L. Fox, and “Aeronomy” (chapter for the 
Atomic and Molecular Physics Reference Book, AIP Press, 1996, edited by Gordon Drake) 
by J. L. Fox. 

In addition, two talks based on work supported by this grant were given at annu- 
al meetings of the Division of Planetary Sciences of the American Astronomical Society: 
“Heating efficiencies in the Titan thermosphere” (J. L. Fox and R. V. Yelle) Palo Alto, 
November, 1991; “A new model of the ionosphere of Titan” (J. L. Fox and R. V. \elle), 
Kona, October, 1995. Invited talks were also presented by J. L. Fox at the Second Interna- 
tion Symposium on Dissociative Recombination (“Dissociative recombination in planetary 
ionospheres”) in St. Jacut. May, 1992, and at the Third Internation Symposium on Dis- 
sociative Recombination (“Hydrocarbon ions in the ionospheres of Jupiter and Titan”) in 
Ein Gedi, May, 1995. The subject matter of these talks included discussion of the Titan 
ionosphere based on work supported by this grant. 
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Heating Efficiencies in the Thermosphere of Titan 
J. L. Fox 

State University of New York at Stony Brook 

R. V. Yelle 
Boston University 


Abstract 

We have constructed a coupled ion and neutral chemical model of the ionosphere 
and thermosphere of Titan. Along with density profiles of ions and minor neutrals, we 
have computed the heating rates and heating efficiencies for the neutrals and chemical 
heating rates and efficiencies for the ions. We find that the neutral heating efficiency in 
our standard model varies from about 30% near 800 km to 22% near 2000 km. The most 
important heating processes are neutral-neutral reactions and photodissociation. The ion 
chemical heating rates maximize at about 10 eV cm^V 1 , and the corresponding ion 
heating efficiencies peak near 1000 km at about 0.6%. 
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Introduction 

The heating efficiency is the fraction of solar energy absorbed that is deposited locally as 
heat. The EUV (~100-1000 A) and FUV (~1000-2000A) regions of the solar spectrum 
are absorbed in the thermospheres of planets. Solar energy is transformed to heat in pho- 
todissociation and photoelectron impact dissociation of molecules, in exothermic reactions, 
including ion-molecule reactions, neutral-neutral reactions, and dissociative recombination 
of ions with electrons. Neutral heating efficiencies are important in determining the neu- 
tral temperature structure and thus the chemical structure of the thermosphere. Although 
the main source of ion heating in non-auroral planetary ionospheres during the daytime is 
Coulomb collisions with photoelectrons, ion-molecule reactions are a source of heat that 
can be important in the region of the ion peak (Rohrbaugh et al., 1979; Cravens et al., 
1980). 

We have modeled the thermosphere and ionosphere of Titan, with a view toward 
determining the ionospheric structure, altitude profiles of minor neutral species, and the 
EUV heating efficiencies. We have adopted the F79050N solar fluxes of Hinteregger ipri- 
vate communication; see also Torr et al.. 1979], which is appropriate to a period of high 
solar activity, such as that which prevailed at the time of the Voyager 1 encounter with 
Saturn and Titan (1981). The solar zenith angle in our model is 60°. Photoionization 
and photodissociation are included for No, CH4, C2H6, C2H4, C2H2, and H2, as well as 
photoionization of H, N and C. In addition, photoionization of CH 3 by Lyman alpha, and 
photodissociation of HCN are also included. Electron impact excitation and ionization 
are included for N 2 , CH 4 , H 2 , H. C and N. The cross sections for photoionization and 
dissociative photoionization of N 3 from threshold to 116 A were taken from Samson et al. 
i 1987] . The N 2 electron impact cross sections are the same as those adopted by Fox and 
Victor 1 1988 1 in their study of electron energy deposition in N2 gas. The cross sections for 
photoionization of hydrocarbons, and electron impact on H2, H and CH 4 are the same as 
those compiled by Kim and Fox 1990 . The cross sections for photodissociation of HCN 
were taken from Nuth and Glicker 1982] and Lee 19801. The cross sections for photoion- 
ization of C and N were taken from LeDourneuf et al. [1976], Cantu et al., 1981 . and 
Daltabuit and Cox [1972]. The electron impact ionization cross sections for C and N were 
adopted from the measurements of Brook et al. 1978]. The ionization potential of CH 3 is 
9.84 eV, and photoionization by solar Lyman alpha is possible, but the cross section has 
not been computed or measured. We assume, as did Strobel [1975], that the cross section 
is equal to that of atomic carbon at threshold. 1.2 x 10~ 1 ' cm 2 . 
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Sources of Heating 

Photodissociation . In photodissociation, the amount of energy that goes into kinetic energy 
is the difference between the energy of the photon and the dissociation energy (which may 
include some internal excitation of the fragments). In the absence of information about 
the electronic states of the fragments, they are assumed to be produced in their ground 
states. When there are molecular products, some of the energy may be taken up as 
vibrational energy of the fragments. The energy that appears as vibrational excitation in 
photodissociation has been found to be small [cf. Fox, 1988 and references therein!; this 
would be expected to be especially true when the fragments are H atoms. We assume here 
that 25% of the excess energy appears as vibrational energy of molecular products, and 
the remainder is converted to kinetic energy, as we have in previous calculations of heating 
efficiencies je.g. Fox, 1988). A half- collision model (Holdy et al., 1970) suggests, however, 
that when one of the fragments is very light, such as a hydrogen atom or molecule, the 
energy that appears as vibrational excitation of the other fragment may be very small. 
It is possible that only 10-15% of the energy is taken up as vibrational excitation in 
photodissociation. 

The N 2 dissociation threshold is large, about 9.76 eV, and because most of the dissoci- 
ation occurs through discrete excitation into predissociating states in the singlet manifold, 
the effective dissociation threshold is larger. Consequently, N 2 does not absorb longward of 
about 985 A. Although N 2 is the major species in the atmosphere, the peak heating rate due 
to photodissociation of N 2 is only 5.7 eVcm -3 s _1 at about 1060 km. near the ionospheric 
peak. By contrast. CH4 is a strong absorber, and photodissociation of CH 4 is the most 
important heating process of this type. The peak heating rate is about 1.4 x 10 3 e\ cm 3 s 
at 770 km. Although at this altitude infrared heating may be more important than solar 
ultraviolet heating, the dominance of heating due to methane photodissociation extends 
to higher altitudes. At 1060 km. the heating rate due to methane photodissociation is 
1.23 x 10 2 eVcm -3 s -1 . Even HCN photodissociation is more important than that of N 2 , 
with a heating rate of about 8 e\ cm -3 s _1 at 1060 km. 

Chemical reactions. The energy released in exothermic reactions can appear as ki- 
netic energy of the products, thus resulting in heating. The partitioning of kinetic energy 
between the products in any reaction is determined by conservation of momentum and 
energy. For the case in which there are only two products, the kinetic energies of species 
1 and 2, £1 and £ 2 , respectively, are given by 

it = tn 2 /(m 1 - m 2 )e 

and 

£ 2 — ttl \ / ( rn ] — m 2 )£ 

where m ; is the mass of species j and £ is the total kinetic energy released in the reaction. 

Dissociative recombination reactions in particular are often very exothermic and con- 
tribute significantly to neutral heating. On Mars and Venus, dissociative recombination 
of the major ion, 0,", is the most important chemical heat source for the neutrals Fox 
and Dalgarno, 1979, 1981; Fox, 19881. On the Earth, dissociative recombination of NO + 
is the most important heat source among reactions of species in their ground states [Torr 
et al., 1980). On Titan, we find that dissociative recombination is not as important as it 
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is on the terrestrial planets, partly because the ions are “processed” more, that is. they 
have undergone more transformations before they recombine. More important, however, is 
the presence of methane, which is a strong absorber, increasing both the heating rate due 
to photodissociation and the production rate of neutral fragments that can participate in 
neutral-neutral reactions compared to the terrestrial planets. 

Dissociative Recombination. In dissociative recombination of polyatomic molecular 
ions, as in photodissociation, some uncertainty exists about how much of the exothermicitv 
appears as vibrational excitation of the molecular products. And just as for photodisso- 
ciation, a half- collision model suggests that little energy is expected to be deposited as 
vibrational excitation. We have assumed that 25% of the energy released appears as vi- 
brational energy of molecular fragments; that value may be an overestimate where one of 
the products is an H atom. For example, an important dissociative recombination reaction 
for neutral heating is that of HCNH + . There are four exothermic product channels for 
this reaction: 

HCNH+ + e -» HCN + H 
-» HNC + H 
— ► CH 2 + N 

- h 2 ^cn. 


The only information that is available about the product channels is the measurement of 
Adams et al. 1991 i , that indicated that 0.63 H atoms are produced per recombination. 
This implies that the sum of the branching ratios for the first two channels is 63% and 
the sum of the second two branching ratios is 37%. We somewhat arbitrarily assume that 
production of HCN, the more stable isomer, is the dominant channel of the first two. and 
that the third channel is more important than the fourth. 

Electron impact heating. In electron impact dissociation some of the energy lost by the 
electron may appear as kinetic energy of the neutrals produced. In general, the fraction 
of energy that appears as kinetic energy is not fixed for a given process and it is therefore 
not possible to determine the heating due to electron impact dissociation exactly. An 
average value may be determined from time-of-flight studies, or derived from electron 
energy loss spectroscopy, but measurements are not available for all of the electron impact 
dissociations. Fortunately, Prokop and Zipf [1982] have measured the average kinetic 
energy of the N atoms that are released in electron impact dissociation and dissociative 
ionization of N 2 , and report values of 0.45 eV and 3.0 eV, respectively. Since this is the 
most important species in the Titan atmosphere, it accounts for most of the photoelectron 
energy loss at the ion peak, near 1000 km. In computing the heating due to electron 
impact dissociation, we have adopted the measured values for N 2 , and we have assumed a 
value of 1.0 eV per dissociation or dissociative ionization of H 2 and CH 4 . 

Ion-molecule reactions. In ion-molecule and neutral-neutral reactions, the exother- 
micity can appear as internal energy of the products, as well as kinetic energy. In general, 
little information is available about the partitioning of energy among the various modes. 
For our standard calculation, we have assumed that 60% of the energy appears as vibra- 
tional energy of the products. We also then compute lower and upper limits to the heating 
efficiency using the values 80 and 40%, respectively. As mentioned above, for ion- molecules 
reactions, the partitioning of kinetic energy between the the neutral and the ion product 
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is determined by momentum conservation. The most important ion-molecule reactions for 
thermospheric heating are 

NT - CH 4 -> CH; - N 2 -r H 

and 

ch~ -ch 4 -» c 2 e; ~h 2 . 

The most important neutral-neutral reactions are: 

CH 2 t CH 2 -*• C 2 H 2 tH 2 


and 


XHrN^N 2 -H. 


Results 

The neutral heating rates as a function of altitude are presented in Figure 1 and resulting 
heating efficiencies are shown in Figure 2. The standard values vary from about 0.30 at 
800 km to 0.22 at 2000 km. The lower and upper limits are smaller and larger by about 
0.05 near 800 km but converge to nearly the standard value at 2000 km. These heating 
efficiencies are very different from those appropriate to Earth, which are in excess of 0.50. 
and are more like those for \ enus and N'lars, which are about 0.20 - 0.24 It is apparent that 
in N 2 atmospheres, it is the minor species that dominate the absorption and heating: 0 2 
for the terrestrial case, and methane for Titan. 

The ion chemical heating rates are shown in Figure 3 for the standard, upper and 
lower limit models and and the corresponding heating efficiencies are shown in Figure 4. 
The chemical heating rates are much smaller than the neutral heating rates, and attain 
maximum values on the order of 10 e\ cm 3 s ^ near the ion peak. The heating efficiencies 
maximize near the ion peak at values of 0.003, 0.006. and 0.009 for the lower limit, stan- 
dard, and upper limit models. We can compare the ion chemical heating rates to those 
computed by Roboz and Nagy (1994). Their "solar-onlv” heating rates maximize near 
the ion peak with values that are much less than ours, about 0.3 eVcm 3 s 1 . The source 
of the difference is uncertain, but probably arises from our use of a more complete set of 
reactions. 
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Figure 1. Altitude profiles of the heating rates of the most important processes in the 
Titan thermosphere. 


Figure 2. Altitude profiles for the heating efficiencies for solar radiation in the 14-2000A 
interval. The standard, upper limit and lower limit profiles are illustrated. 


Figure 3. Altitude profiles of the ion chemical heating rates in the Titan thermosphere 
for the standard (solid curve), lower limit (short-dashed) and upper limit (long-dashed) 
models. 


Figure 4. Altitude profiles for the ion heating efficiencies for solar radiation in the 14- 
2000A interval. The standard (solid curve), upper limit (short-dashed) and lower limit 
(long-dashed) models are illustrated. 
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Abstract 

We have constructed a new model of the ionosphere of Titan that includes 69 species and 
626 reactions. Although NJ is the major ion produced in the ionosphere, the ionization 
flows to ions whose parent neutrals have lower ionization potentials and to species with 
large proton affinities. In contrast to other models, which have predicted that HCNH 
should be the major ion, we suggest that the major ions at and below the ion peak near 
1050 km are hydrocarbon ions and H, C, and N-containing ions. Our predicted peak 
electron density for a solar zenith angle of 60° is 6.3 x 10 3 cm -3 at an altitude of 1050 km. 

Introduction 

In addition to N 2 and CH 4 , the thermosphere of Titan contains significant densities of H 2 , 
H, and non-methane hydrocarbons such as acetylene, ethane, and ethylene (e.g., Hunten 
et al., 1984, Strobel et al., 1992), and appears to be well-mixed to very high altitudes. The 
atmosphere may also contain as much as 14% Ar. The homopause in the model of Strobel 
is near 1050 km, where the eddy diffusion coefficient is about 1 x 10 9 cm 2 s -1 . The exobase 
is near 1500 km. and the density profiles of H and H 2 are altered appreciably by thermal 
escape at the top of the atmosphere. 

Information about Titan’s ionosphere is limited to the radio occultation measurements 
made by the Voyager 1 spacecraft. Lindal et al. (1983) reported only upper limits of 
(3 — 5) x 10 3 cm -3 near the terminators as the spacecraft entered and exited the occultation 
region. Recently. Bird et al. (1995) have reanalyzed the Voyager radio occultation data 
and reported a possible detection with a peak electron density of 2700 cm -3 at an altitude 
of 1190 km near the evening terminator. 

The ionosphere of Titan has been predicted to arise from both solar photoionization 
and from the interaction of the neutral atmosphere with energetic electrons from Saturn's 
magnetosphere, in which the orbit of Titan is embedded during part of its orbit (e.g.. 
Strobel and Shemanskv, 1983: Atreya, 1986). Strobel et al. (1991) found, however, that the 
N + 1085 A and X 2 Lyman-Birge-Hopfield band intensities that were measured by Vovager 
1 spacecraft during its Titan flyby in 1980 could be reproduced with the solar source alone. 
Cravens and co-workers [Keller et al., 1992; Gan et al., 1992] constructed a model of the 
ionosphere of Titan that included the interaction of energetic magnetospheric electrons 
with the thermosphere and found it to be a minor source. We model here the ionosphere of 
Titan produced by photoionization and photoelectron-impact ionization alone. Although 
we ignore here the interaction of the thermosphere with Saturn’s magnetosphere, we do 
not mean to imply that magnetospheric electrons are unimportant all the time or for every 
process. In all probability, however, they are unimportant for ion production near and for 
a substantial distance above the ion peak. 
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The Model 

The background model neutral thermosphere that we have adopted here is shown in Fig. 
1. The density profiles of N2 and CH4 are taken from Strobel et al. (1992), and mixing 
ratios of C 2 H 2 , C 2 H 4 , C 2 H6, and C 4 H 2 are taken from Yung et al. (1984). The mixing 
ratio of methane has been estimated as 2-10% in the Titan atmosphere (Hunten et al.. 
1984), but we have adopted a mixing ratio of 3% at the lower boundary, consistent with 
the preferred model of Strobel et al. (1992). The CO mixing ratio is adopted from the 
ground-based measurements of Gurwell and Muhleman (1995). The temperature profile is 
taken from Strobel et al. (1992), and is nearly isothermal at the exospheric temperature 
of approximately 175K above about 1000 km. The electron temperatures are taken from 
Gan et al. (1992). 

In our model, the density profiles of six species were fixed, and we compute the densi- 
ties of 34 neutral species and 35 ions, which are enumerated in Table 1. We have included 
photoabsorption and photoionization of N 2 , CH 4 , C 2 H 2 , C 2 H 4 , C 2 H6, N, C, H and H 2 , pho- 
todissociation of HCN, and photoionization of CH3 by Lyman alpha. The solar spectrum 
used is the high solar activity F79050N spectrum of Hinteregger (private communication: 
see also Torr et al.. 1979), and the solar zenith angle is 60°. Photoelectron-impact ion- 
ization was included for X 2 . CH 4 . H. H 2 , N and C. The most important production rate 
profiles are shown in Fig. 2. The major ion produced is NT from 600 to 1800 km. and 
CH'J’ above that altitude. Other important ions produced include N~. CH 3 , C 2 H., . and 
CoHT. 

The ion density profiles were computed including 626 chemical reactions. The ion- 
molecule, neutral-neutral and dissociative recombination reactions in the model are listed in 
Table 2 along with their assumed rate coefficients. Most of the rate coefficients for the ion- 
molecule reactions were taken from the compilations of Anicich and Huntress (1986) and 
Anicich (1993). For neutral-neutral reactions, rate coefficients were taken from the model of 
Yung et al. (1984) and from models of the interstellar medium (e.g., Millar et al., 1988). For 
many of the most important dissociative recombination coefficients, the rates are unknown: 
where no measurements are available, values of (3 — 3.5) x 1CP ' (300 /T f ) U o cm 3 s 1 were 
assumed. Even where dissociative recombination coefficients have been measured, in many 
cases there is little information about the product yields, for which in general there are 
many possibilities, especially for polyatomic ions. Because the ion and neutral chemistries 
are coupled, (and because the exothermicity and production of fast atoms depend on the 
channel by which the recombination proceeds), measurements of the product yields are 
important. Here we have assumed that the available channels are populated with equal 
probability. Eddy and molecular diffusion for neutrals, and ambipolar diffusion of ions 
were included in our model, for which the lower and upper boundaries were at 600 km and 
3000 km, respectively. 
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The predicted densities of the major ions are shown in Fig. 3. The peak electron 
density is about 6.3 X 10 3 cm -3 near 1050 km. This is larger than either the Voyager 
upper limits and the possible detection of Bird et al. (1995), but those measurements are 
for the terminator region. We would expect the densities to be substantially larger than 
ours, which are for the davside at 60° solar zenith angle. 

Although is the major ion produced up to about 1800 km, below about 1900 km 
it is lost by reaction with methane 

N ' - CH 4 — CH+ + N 2 + H 2 
-> CH 3 + N 2 + H, 

and above that altitude by the reaction: 

-H 2 -4 N 2 H+ + H. (2) 

Reaction with H and charge transfer to acetylene are minor loss processes in the lower iono- 
sphere. Dissociative recombination becomes the second most important loss mechanism 
above about 2500 km. Altitude profiles for the loss processes are shown in Fig. 4. 

The major ion at high altitudes is N 2 H + . The major production and loss mechanisms 
for N 2 H + are shown in Fig. 5. X 2 H~ is produced mainly in the reaction of NT with H 2 
(reaction 2) below 1300 km and above 1700 km. Over a limited altitude range, about 
1300-1700 km. the most important source is 

HT - N 2 NolT - H (3) 


Other minor sources include the reactions of H 3 with N 2 and NH with N 2 . The most 
important loss mechanism is dissociative recombination 


X 2 H r — e — ♦ products 

(4) 

above 2300 km and proton transfer to CH 4 : 


N 2 H* -CH 4 - CHt + N 2 

(5) 

below that altitude. Reaction with C 2 H 2 : 


n 2 h* + c 2 h 2 -> c 2 h+ + n 2 

(6) 


is also significant near the ion peak. 

In the altitude region from 1300-1800 km, the major ion is C 2 H7. At the peak near 
1050 km the dominant ion is the pseudo-ion C X H+ . The density of C 4 H^ exceeds that of 


(la) 
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CzHjJ" near its peak at about 750 km. Most other models have predicted that HCNH~ 
would be the major ion (Ip, 1990; Keller et ah, 1992; Keller and Cravens. 1994). The 
major production and loss mechanisms for HCNH + in our model are shown in Fig. 6. 
Above about 1050 km, HCNH + is produced in the reaction; 

N + — CH 4 — HCNH^ + H 2 . (7) 

Below the ion peak the most important source is proton transfer from C2H 5 to HCN: 

C 2 H7 -f HCN — * HCNH + 4- C 2 H 4 . (8) 

Minor sources include the reactions of HCN + ■with CH 4 at high altitudes and proton 
transfer of C2H3" to HCN at low altitudes. Below the ion peak, HCNH~ is lost mainly by 
proton transfer to C 4 H 2 : 

HCNH~ - C 4 H 2 -» C 4 H+ + HCN. (10) 

The proton transfer reactions with HC3N, and with NH3 are also significant below about 
1000 km. Loss is by dissociative recombination above that altitude. 

The major production reactions for C r Hj are reactions of C 4 - and C 2 - hydrocarbon 
ions with acetylene. C 2 - hydrocarbon ions with methane and a few reactions of Cj- or C 2 - 
hydrocarbon ions with ethylene and ethane. The most important production reactions are 


plotted in Fig. 7 and are listed below: 

CH^ -C 2 H 2 - C3H3- ~H 2 (10) 

c 2 h; -c 2 h 2 - c 3 h* tCH 3 (11) 

C 2 H; -C 2 H 2 — C 4 H+ + H (12) 

C 2 HT - C 2 H 2 - C3H+ + CH 4 (13) 

C 2 Hr - C 2 H 2 - C 4 Ht + H 2 (14) 

C 2 H7 - CH 4 - C3 H+ +H (15) 

C 2 H* -CH 4 - C3H+ + H 2 (16) 

C 2 H~ -C 2 H 4 - C 3 H^ +CH 3 (17) 

C 2 H7 A C 2 H 4 — > C 3 Ht + CH 4 (18) 

CHJ + C 2 H6 - CaHf +H 2 (19) 


Although C r Hy clearly does not represent a single ion, and thus is not the major ion 
in that sense, the computed peak density of HCNH + is almost an order of magnitude less 
than the total electron density. 
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Photoelectron fluxes 

A secondary but important source of ions near and below the ion peak is photoelectron 
impact ionization. The predicted primary and steady-state fluxes from of photoelectrons 0 
to 200 eV at 0.5 eV resolution near 1000 km are shown in Fig. 8. The top spectrum is the 
initial spectrum and the bottom is the steady-state spectrum averaged over three intervals 
(1.5 eV). The spectrum displays a sharp drop-off above 70 eV, where photoelectrons pro- 
duced by the the last of the strong solar lines appear. The spectrum appears similar to 
that of Gan et al. (1992) for their solar-only case, but the sharp bite-out from about 20 to 
40 eV does not appear in our spectrum. The reason is unknown, but will be investigated. 

Discussion 

Solar ionizing photons are absorbed in the thermospheres of planets and satellites, and the 
peak ionization rate is usually found at a column density of about 10 1 ‘ cm 2 . This altitude 
may be near or above the homopause, where the transition from transport by large scale 
mixing to that by molecular diffusion takes place. Although the density of CH 4 eventually 
exceeds that of N 2 at high altitudes on Titan, the ionization peak appears roughly at the 
homopause, which is near 1050 km in the Strobel (1992) model. There the major neutral 
species is N 2 and the mixing ratio of CH 4 is about 3%. Although hydrocarbon and related 
ions clearly comprise only a small fraction of the ions produced by photoionization and 
electron-impact ionization on Titan, our model shows that such ions may be abundant 
near and below the major ion peak. 

The major ions produced in any ionosphere will be transformed by ion-neutral re- 
actions if they are formed in the presence of sufficient densities of neutrals with which 
they can react. In reducing environments, such as the atmospheres of the outer planets 
and their satellites, ionization tends to flow from species whose parent neutrals have small 
proton affinities to species whose parent neutrals have large proton affinities. On the giant 
planets, the terminal ions will be large protonated hydrocarbons (e.g., Kim and Fox. 1991; 
1994); on Titan, the major species will be hydrocarbon ions, protonated nitriles and other 
ion containing C. H and N. In most cases, the available information about ion chemistry 
is insufficient to identify the terminal ions in the bottomside ionospheres, but since larger 
molecules tend to have larger proton affinities, the average size of the ion probably increas- 
es with decreasing altitude. In order to improve our understanding of these ionospheres, 
rate coefficients and branching ratios for reactions of hydrocarbons and hydrocarbon ions, 
and neutral and ionic species containing C, H and N are needed, including those for disso- 
ciative recombination reactions. Some representative proton affinities are shown in Table 
3. 

Clearly, a major uncertainty in the ionosphere model is the accuracy of the neutral 
model, including the density profiles computed by Yung et al. (1984). We note here that 
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the crucial C 4 H 2 density was underestimated in the middle atmosphere in the Tung et 
al. model compared to measured values. It seems that, if anything, the actual values 
will be larger than those we have assumed. Moreover, the conclusion that HCNH does 
not dominate the ion density profile depends only on the presence (with a large enough 
density) of some species that has a higher proton affinity than HNC or HCN, or of some 
other chemical pathway to the protonated species. With a fairly large thermospheric 
mixing ratio of methane, and the presence of solar ultraviolet photons, that is all but 
certain. 
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Figure Captions 


Fig. 1. Neutral density profiles in the thermosphere of Titan that were adopted in these 
calculations. The densities of X 2 . CH 4 , C 2 H 2 , C 2 H 4 , C 2 H 6 , and CO were assumed to be 
fixed. The density profiles of H 2 and H were computed self-consistently in the model. 

Fig. 2. The computed production rate profiles for the most important ions produced in 
the Titan ionosphere. 

Fig. 3. Computed steady-state densities of the most important ions produced. The solid 
curves are labeled by the ion produced dashed and the dashed curve is the electron 
density profile. C r H+ is a pseudo-ion that represents all hydrocarbon other than those 
included explicitly. 

Fig. 4. Altitude profiles of the major loss processes for N^" in the Titan ionosphere. The 
curves are labeled by the species with which the ions react. 

Fig. 5. Altitude profiles for the major sources (top) and sinks (bottom) of N 2 H + in the 
Titan ionosphere. 

Fig. 6. Altitude profiles for the major sources (top) and sinks (bottom) of HCNH + in the 
Titan ionosphere. 

Fig. 7. Altitude profiles of the major sources of C r HJ in the Titan ionosphere. 

Fig. 8. Computed primary (top curve) and steady-state (bottom curve) photoelectron 
fluxes at 1050 km. The bin size for the primary fluxes is 0.5 eV. The steady-state fluxes 
have been averaged over three bins. The peak at 200 eV is an artifact that results 
from the treatment of photoelectrons with energies larger than 200 e\ . which have been 
divided into one or more 200 eV electrons and an electron with the remainder of the 
energy. 
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1.350 x 10" 9 
1.350 x 10 -9 

8.800 x 10-'" 

6.160 x 10- |M 

1.100 x 10- 9 

1.800 x 10" 9 
3.680 x 10-"' 

1.230 x 10-" 

5.160 x 10" 10 

5.880 x 10-"' 

3.750 x lO" 10 
7.950 x 10- 10 

3.300 x lO - 10 
2.630 x lO" 10 

1.310 x lO" 10 
8.760 x 10-" 

1.460 x 10-" 

7.880 x 10“ 10 

7.300 x lO”" 

1.310 x lO -10 
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C,HT+ N — CH+4- HCN -f 
C 2 H7 + N — C,N+4 H 2 + 

C 2 H j + N — C X H V N+- H - 
C 2 H7 + NH 3 — NHJ+ C 2 H 2 + 

C 2 Ht + NHs— NH|- C 2 H+ 

C 2 Ht + HCN — H 2 CN+ + C 2 H- 
C 2 H++ HCN — C X H,,N+- H - 
C 2 H++ HCN — C x H, y N+- 4- 
C 2 Ht+ HC 3 N— » C X H+ + HCN - 
C 2 H.t+ HC 3 N^ C.H,N++ J- 
C 2 Ht+ H 2 — C 2 HJt H + 

C 2 H7 + CH 3 NH 0 — • C x H y N+- C 2 H- 
C 2 H++ CH 3 NH 2 ^ C r H y Nt 4 - C 2 H 2 - 
C 2 H++ CH 3 NH 2 — C,H f N++ C 2 H 3 ^ 
C 2 Ht+ CH 3 CN— C x H y N+ + C 2 H- 
C 2 H:4 CH 3 CN— C x H+ + CN - 
C 2 Ht+ CH 3 CN— C x H+ + HCN - 
C 2 Ht + CH 3 CN— C x H y N++ 4- 
C 2 Ht+ H 2 ^ C 2 H|+ - 


c 2 n++ 

HCN - 

C X H V N++ - 

C,H+ 

4 

NH — C X H„N++ H - 

C 2 HJ 

4 

CH 4 ~ ► 

: x h+- 

H 2 -r 

C,H : | 

4 

c 2 h 2 — 

C ^ H v- 

H 2 -r 

c 2 h+ 

4 

c 2 h 2 — 

C 2 H;7- 

C 2 H 2 - 

CjH r | 

4 

C 2 H 4 - 

CiH7 — 

c 2 h 2 - 

c,h'| 

4 

c 2 h 6 - 

C 2 H 7 - 

c 2 h 4 * 

CjH.{ 

4- 

c 2 h«- 

c x H r ;- 

ch 4 - 

CjH+ 

4 

c 2 h 6 - 

c x h.;- 

h 2 - 

c 2 h+ 

4 

HCmN — 

■ C X H V > 

IT- C 2 K 2 - 

c,h; 

4 

c 2 n 2 -> 

C r H„N 

7 - CjH 2 — 

Cjh': 

- 

c 2 n 2 -> 

C x H y N 



c-,h.t 

4 . 

HCN — 

H>CN* 

c,h 2 - 

C.HJ 

4 

c 2 n 2 -> 

C X H V N 

7 ~ C 2 Hj- 


C.H++ H — C 2 HT-r H 2 — 

C,H-|+ N — C x H y N7- H 2 - 

c.h; - c 2 h 2 — c x h*- ch 3 - 
c 2 h++ c 2 h 2 — c x h ! ;- h - 

C.H++ C 2 H 4 — C x H~- CH,- 

c 2 h++ c 2 h 4 — c x h;- h - 
c 2 h++ c 2 h 6 — C x H j " 4 ch 4 - 
C 2 H++ C 2 H 6 - C x H+4- ch 3 - 

c 2 h++ nh 3 — nh 3 + c 2 h 4 ^ 

C 2 H++ NH 3 — NHJ-r C 2 H 3 + 
C 2 H+ + HCN — C x H y Nt~ 4- 
C 2 H+-f HC 3 N— C,H,N++ CjHs- 
C 2 H++ H — C 2 Hj-T h 2 + 

c 2 h++ ch 4 ^ c x h+- h 2T 
c 2 h++ C 2 H 2 - C x H+~ CH 4 ^ 
C-.H++ C 2 H 2 — C x H + + H 2 + 
c 2 h++ c 2 h 4 -> C x H^4 CH 4 + 

C 2 H++ C 2 H 6 - C x H+~ H 2 -r 


2.500 x 10~ il 

7.500 x 10- 11 

1.500 x lO" 10 
1.240 x 10~ 9 
1.860 x 10“ 9 
6.720 x 10~ u 
3.780 x 10' 11 

5.000 x 10 _28 (T n /300) o,oo exp(— 0.00/T n ) [M: 

1.480 x 10- 9 

2.220 x 10~ 9 

1.000 x 10 -11 

1.300 x 10- 9 
7.560 x 10-'° 

6.480 x 10- 10 
8.360 x 10~ l ° 

1.060 x 10“ 9 

1.060 x 10- 9 

5.000 x 10-- 8 (T„/300) o oo exp(-0.00/T n ) [Ml 

1.200 x 10- 27 (T„/300)°- oo exp(-0.00/T n ) [Ml 

4.200 x 10" 10 

6.500 x lO" 10 

2.000 x 10” 10 
2.160 x 10" 10 
5.040 x 10" 

9.300 x 10“ 111 
2.910 x 10 _i " 

2.480 x 10" 10 

8.060 x 10" u 
3.800 x lO" 9 

5.500 x 10-'" 

1.000 x 10- JS (T n , 300V’°"expt— 0.00 T n ) >1] 
1.600 x lO -9 

4.000 x 10- 1,1 

1.000 x 10 ""’ 

2.200 x 10" 1,1 
6.730 x 10- 10 
2.370 x 10- 10 
7.190 x lO” 10 
7.110 x lO" 11 
3.710 x lO" 13 
4.930 x 10 -12 
1.240 x 10“ 10 
1.940 x 10 -9 

1.000 x 10 -28 (T„ /300) lU,o exp(— 0.00/T„) [Ml 
1.100 x 10” 9 

3.000 x 10” 10 

1.000 x lO" u 
6.840 x 10- u 

1.220 x lO” 10 
3.550 x lO” 10 
3.270 x 10" u 
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C 2 Hj+ C 2 H 6 - C z H + -r CH 4 - 
c 2 h++ nh 3 — NH+ + c 2 h 4 - 
C 2 H++ CH 3 NH 2 — C x H,N+- H 2 — 
C.Ht+ HC 3 N— * C r H„N+- C 2 H 4 - 
c 2 r ++ hcn — H 2 CN+- C 2 H 4 - 
C 2 H^+ CH 3 CN— C x H y N+- C 2 H 4 - 
C 2 Ht+ C 2 N 2 — C x H y NT + C 2 H 4 - 
C .H++ C 2 N 2 - C r H,N+-i- - 
Ci Hg + C 2 H 2 — C 2 Ht+ C 2 H 3 - 
C 2 H++ C 2 H 2 - C x H+ + CH 3 - 
C 2 H++ C 2 H 2 - C x H+ + H - 

c 2 h++ c 2 h 4 — c 2 h|+ c 2 h 6 - 
c 2 h++ c 2 h 6 - c x h++ ch 4 - 
c 2 h++ c 2 h 6 - c x h++ ch 3 - 

C 2 H++ NH 3 - NH++ C 2 H 6 - 

c 2 h++ nh 3 - NH++ c 2 h 5 - 
c 2 h++ h — c 2 h^+ h 2 + 
c 2 h++ CH 4 - c 2 h++ + 
c 2 h++ c 2 h 2 — c 2 hj+ c 2 h 6 - 
c 2 h++ c 2 h 4 - c 2 h++ c 2 h 6 - 

C 2 H++ NH 3 — NH+ + C 2 H 6 - 
C.H++ HCN — H 2 CN+- C 2 H„- 
C 2 H|+ HCN — Cj-HyN?"— CH 4 - 
N++ H 2 — NH+4- H + 

n++ ch 4 — ch;- nh - 
n + + ch 4 — ch;- n - 
N+J- ch 4 — h 2 cn + - h>- 


N- + 

ch 4 ^ 

HCN T — Ho-r H 

N r + 

c 

2 h 4 - 

C-jHT-f- NH 2 - 

N + + 

c 

-h 4 - 

C 2 H~h- NH - 

N+ + 

c 

-h 4 - 

c 2 h;+ n - 

N r + 

c 

»h 4 — 

HCN*+ CH;,- 

N+4- 

c 

2 H 4 — + 

H>CN+- ch 2 - 

N+4- 

c 

2 h 4 — 

C x H,,N+~ H 2 - 


N + + NH;,— NH 3 -r N - 
N + + NH;,— NHT - NH - 
N++ NH 3 — N 2 H + - H 2 - 
N"*"+ CH;,NH 2 — C x H„N+- N - 
N + + CH;,NH 2 — C X H V N+- NH - 
N + + CH;,NH 2 — C x H y N+- NH 2 - 
N+ + CH 3 NH 2 — H 2 CN + + NH 3 - 
N+ + CH 3 NH 2 — CH 3 + N 2 — H 2 
N + 4- C 2 N 2 — C x H y Nj-+ N - 
N + + C 2 N 2 — C 2 N + + N 2 + 

N+ + HC 3 N- C x H + + N 2 + 

N + + HC 3 N— C x H y N+- N - 
N++ CO — CO + i N 4- 
N + 4- CO — NO+- C + 

N + + NH — NH + 4- N 4- 
N + + NH — Nt+ H 4- 
N + + HCN — HCN++ N + 

NH+4- HCN — H 2 CN + + N - 


5.320 x 10~ 12 
2.180 x 10“ 9 
1.870 x 10“ 9 
3.300 x 10- 9 
2.700 x 10~ 9 
3.800 x 10~ 9 
8.400 x 10 _n 

1.000 x 10- 28 (r n /300) ooo exp(-0.00/T n ) [Ml 

2.220 x 10“ 10 

8.190 x lO" 10 

1.290 x IO -10 

1.150 x 10- 9 

7.980 x 10- 12 
1.100 x 10- 11 
6.240 x 10- 10 
1.610 x 10- 9 

1.000 x 10 -10 

2.000 x 10- 8 (T„/300) o - ooo exp(-84.00/T n ) 
2.000 x 10" 9 

2.000 x 10- 9 
2.000 x 10" 9 

1.980 x 10" 9 
2.200 x 10“ 10 

4.160 x 10-'°exp(-41.90/T n ) 

4.260 x 10- 10 
4.920 x lO” 11 
2.710 x lO" 

7.380 x lO -11 


1.600 

X 

10-"’ 

4.800 

X 

10-1" 

4.000 

X 

10-1" 

2.400 

X 

10-1" 

1.600 

X 

10-"’ 

1.600 

X 

10-1" 

1.070 

X 

10- 9 

1.170 

X 

10-"’ 

1.170 

X 

i 0 -i" 

1.400 

X 

10- 10 

1.400 

X 

10- 9 

1.400 

X 

io-"’ 

2.000 

X 

10- 1" 

1.200 

X 

10-"’ 

4.200 x 

10- 10 

9.800 

X 

10- 10 

2.100 

X 

10- 9 

2.100 

X 

10- 9 

8.250 

X 

10-10 

1.460 

X 

10-'° 

3.700 x lO - 10 

3.700 

X 

lO-io 

1.200 

X 

IO" 9 

1.800 

X 

10“ 9 
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NH+ + N — NT- H + 

NH++ NH — NHT+ N + 

NH+ + CO — NO+- H + 

NH++ CO — HCO^- 1 - N + 

NH+ + H 2 — NHT- H + 

NH+ + H 2 — HT - N - 
NH+ + CH 4 — CHT+ N + 

NH+ + CH 4 — NHT -i- CH 3 + 

NH+ + CH 4 — H 2 CN+ + H 2 + H 
NH+ + C 2 H 4 — C 2 HT+ NH 3 t 
NH++ C 2 H 4 — C 2 Hj+ NH 2 + 

NH+4- C 2 H 4 — C 2 HT+ NH + 

NH++ C 2 H 4 — H 2 CN++ CH-)-f 
NH+ + C 2 H 4 — C r H„N++ CH 2 - 
NH+ + C 2 H 4 — C r H y NT + H 2 - 
NH+ + NH 3 — NHT + NH - 
NH+ + NH 3 — » NHT+ N + 

NH+4- N 2 — N 2 H + + N + 

NH+ 4 - CH 3 NH 2 — H 2 CN+4- NH 3 - H 
NH+4- CH 3 NH 2 — C r H y N+4- NH.vr 
NH+ 4 - CH 3 NH 2 — C r H y Nt + NH 2 - 
NH+4- CH 3 NH 2 — C r H y N+4- NH - 
NH+ 4 - CH 3 NH 2 — C x H,N+ + N - 
NHT-i- H 2 — NHT+ H - 
NHT+ CH 4 — NHT- CHa-r 
NHT - c.h 4 — C 2 H|- nh 2 — 

NH+4- Gin.,— C.Hr- NH - 
NHT- C 2 H.(— » C r H y NT4- CH>— 

NHT - C 2 H 4 — C r H y NT — H - 
NHT- NH-, — ♦ NHT- N H > 2 -r 
NHT-r NH 3 — NHT- NH ~ 

NHT 4 - CH,NH 2 — C r H v Nr- NH,- 
NH.T+ CH 3 NH 0 — C r H y NT- NHi — 
NH.+ - CH 3 NH 2 — C r H y N7- NH - 
NHT+ CH 3 NH 2 — NHT- H3CN ~ 
NHT+ NH — NHT- N + 

NH+4- h 2 — NHT- H +• 

NH+4- CH 4 — NHT- CH 3 + 

NH+4- C 2 H 4 — NHT- C 2 H 3 - 
NH+4- NH 3 — NHT- NH 2 — 

NHT+ CH;,NH,— C r H y NT- NH : <- 
NH+4- CH 3 NH 2 — C r H y NT- NH,- 
NHT+ CH 3 NH 2 — NHT + CH2NH2 - 
NH+4- NH — NHT+ N 4- 
NH+4 CH 3 NH 2 — C x H y NT+ NH 3 - 
NH|+ CH 3 CN- C I H y NT*- NH 3 - 
HCN +4 NH — NHT+ CN 4 
HCN +4 CO — C r H y N+4 CO - 
HCN+4- CO — HCO+4 CN + 
HCN+4- CH 4 — H 2 CN + 4 - CH 3 - 
HCN +4 H 2 — H 2 CN + + H + 

HCN+4 C 2 H 2 — C 2 HT+ HCN j- 
HCN+4 C 2 H 2 — C 2 H+4 CN - 


1.300 x 10“ 9 

1.000 x 10" 9 
5.390 x 10- 10 
4.410 x 10 _1 ° 

1.280 x 10~ 9 
2.250 x 10 _1 ° 

9.600 x 10“ n 
1.920 x 10~ 10 
6.720 x 10- 10 
1.500 x 10“ l ° 

3.750 x 10- l ° 

3.750 x 10" 10 

3.000 x 10~ 10 
1.500 x lO" 10 

1.500 x lO -10 

1.800 x 10“ 9 

6.000 x lO" 10 

6.500 x 10"’° 

4.200 x lO -10 
1.050 x lO" 10 

9.500 x 10“ 10 
4.200 x 10" 10 

4.200 x 10“ 10 

2.700 x 10“ 10 

9.200 x 10“ 10 

4.500 x 10" 10 

3.000 x 10-"' 

4.500 x 10-"’ 

3.000 x lO” 1 " 

6.900 x 10-"" 

1.610 x 10" 9 

4.000 x 10-"’ 

1.100 x 10" 9 

4.000 x 10-"’ 

1.000 x 10-"’ 

7.300 x 10“ 19 

3.100 x 10-‘ 2 exp(— 241.00/TJ 

4.800 x 10“"’ 

1.400 x 10“ 9 

4.700 x lO" 9 

9.000 x 10" 10 

6.300 x 10-"’ 

2.700 x lO - 10 

7.100 x lO -10 

2.000 x 10" 9 

1.600 x 10" 9 

6.500 x lO -10 
3.220 x 10 _1 ° 

1.380 x 10-’° 

1.100 x 10" 9 

8.600 x lO" 10 

9.000 x 10"'° 

3.000 x 10" 10 
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HCN + + C 2 H 2 -^ C r H y N++ H 2 - 

3.000 x 10“‘° 

HCN++ NH 3 -> NHJ - HCN + 

1.680 x 10“ 9 

HCN+ + NH 3 — H 2 CN + J- NH 2 + 

8.400 x 10 “ 10 

HCN+ + NH 3 — NH|- CN + 

1.400 x 10 “ 10 

HCN++ HCN — H 2 CN + + CN + 

1.600 x 10“ 9 

HCN+ + HC 3 N— C,H,N++ CN - 

2.210 x 10“ 9 

HCN++ HC 3 N— C I H V N+- HCN - 

2.390 x 10 “ 9 

H 2 CN++ NH 3 ^ NHJ+ HCN + 

2.400 x 10 “ 9 

H 2 CN++ HC 3 N— C,H,N++ HCN - 

3.400 x 10 “ 9 

H 2 CN++ CH 3 NH 2 — C r H y NT+ HCN - 

1.870 x 10 “ 9 

H 2 CN++ CH 3 NH 2 — C,H r N++ - 

2.310 x 10“‘° 

H 2 CN++ CH 3 CN— C x H y NT+ HCN - 

3.800 x 10 “ 9 

H 2 CN++ C 4 H 2 — C 4 HJ+ HCN - 

1.400 x lO * 9 

Nj + CO — CO + - No'-c 

7.400 x 10“" 

NT + CH 4 -> CH 3 -r N\- H 

1.520 x 10 ~ 9 

NT + CH 4 ^ CHT^ N 2 t H 2 

3.800 x 10“'° 

N++ C 2 H 2 — C 2 HT- N 2 + 

4.300 x 10 “ 10 

N++ H 2 — N 2 H"-" H - 

2.100 x 10“ 9 

Nt + N — N++ N 2 - 

5.000 x 10“ 12 

N+ + NH 3 — NHJ — N 2 t 

2.000 x 10“ 9 

Nt+ HCN HCN“- N 2 + 

3.900 x 10“'° 

N+ + CH 3 NH 2 — C r H v N++ N 2 + 

7.200 x 10 “ 11 

N+ + CH 3 NH 2 — C r H v N++ N 2 + H 

8.760 x 10“'° 

NT + CH 3 NH 2 — CH-7- NH 2 t N 2 

2.520 x 10 “ 10 

NT + HC 3 N— C r H y NT- N 2 - 

2.900 x 10 “ 9 

Nt+ C 2 No— C r H v N7- No — 

8.600 x 10“ 10 

N.T+ CH :) CN— ► C f H y NT+ N 2 -r 

3.150 x 10 “ 10 

NT + CH 3 CN^ C r H.,N7 + N 2 + H 

1.370 x 10 “ 9 

NT + CH 3 CN-o C r H v N_7 - No- Hj 

4.200 x 10“ 10 

NT + NH — NH“- No- 

6.500 x lO - 1,1 

NT + H — H + — No- 

1.200 x 10 “'" 

NT 4 - C — C + - No- 

1.100 x 10“ 10 

N 2 H + + CH 4 — CH7- No- 

8.900 x 10“'" 

NoH+- C 0 H 0 — C.HT- No- 

1.410 x 10 “ 9 

NoH + - CoHfi— * C.Hr- No- 

1.300 x 10 “ 9 

N 2 H++ HCN — H 2 CN + 4- No-r 

3.200 x 10 “ 9 

N 2 H + + CH 3 CN— C r H, ( NT - No- 

4.100 x 10 “ 9 

NoH+J- HC 3 N— C x H„N7-r N 2 - 

4.300 x 10 “ 9 

N 2 H + + C 0 N 0 — C r H y NT+ N 2 -i- 

1.200 x 10“ 9 

N 2 H + + CO — HCO*-’N 2 -r 

8.800 x 10 “'" 

CN++ Ho— HCN*— H + 

5.500 x 10 “ 10 

CN+ + H 2 -> C x H y Nr- H - 

5.500 x 10 “ 10 

CN+ + CH 4 -* CH7- HCN - 

4.500 x 10 “ 10 

CN+ + ch 4 — ch;- CN - 

1.350 x 10- 10 

CN+ + CH 4 -o HCN + - CH 3 -t 

1.350 x 10 “ 10 

CN++ CH 4 — H 2 CN“- CHo- 

9.000 x 10“" 

CN+ 4 - CH 4 ^ C x H y NT+ H 2 + 

9.000 x 10“" 

CN++ C 2 H 2 ^ CoHT- CN - 

1.350 x 10 “ 9 

CN++ C 2 H 2 — C r H*NT+ H + 

1.500 x 10“'° 

cn++ C 0 H 4 -* c 2 h;- CN + 

7.000 x 10 “ 10 

CN+ + C 2 H 4 — HCN^-t- C 2 H 3 + 

2.500 x 10“'° 

CN++ C 2 H 4 ^ C x H„NT+ H 2 + 

5.000 x 10 “ u 

CN++ C 2 H 6 - CoH 3 - HCN + Ho 

2.850 x 10 “ 10 


8 


CN+ + C 2 H 6 — C_,H7- HCN - H 
CN++ C 2 H 6 — C 2 Hr- HCN - 
CN++ NH 3 — NHT- HCN 4- 
CN++ NH 3 — NH j - CN - 
CN++ NH 3 — HCN*— NH 2 - 
CN+ + NH 3 — H 2 CN*- NH - 
CN+4- HCN — HCN*- CN - 
CN+ + HCN — C r H„N++ H - 
CN++ HC.,N- C x H v N+ + CN - 
CN+ + HC : ,N- C r H y N++ HCN - 
CN+ + C 2 N 2 — C r H s N7-f. CN - 
CN+4- C 2 N 2 — C x H„N+4- N 2 t 
CN+ 4- C 2 N 2 — C 2 N*- CNN - 
CN+4- CO — CO*- CN - 
CN+4- NH — NH*- CN - 
C 2 N+4- CH 4 — C 2 H 3 - HCN - 
C 2 N+4- CH 4 — H 2 CN‘ r - t - C 2 H 2 - 
C 2 N+4- CH 4 — C,H,N++ H 2 t 
C 2 N+4- C 2 H 2 — C x H*- HCN - 
C 2 N++ C 2 H 2 — H 2 CN~~ - 
C 2 N++ C 2 H 4 — C.H7- - 
C 2 N+4- C 2 H 4 — C x H.,N+4- C.H 2 - 
C 2 N+4- C 2 H 4 — C x H7- HCN - H 
c 2 n+- C 2 H 4 — C X H.;- HCN - 
c 2 n*- C 2 H 4 — C r H y N'r— C 2 H 4 - 
C 2 N+4- NH 3 — H 2 CN*- HCN - 
C 2 N+4- NH.,— N 2 H*- C,H 2 - 
C-.N+4- CH;,CN — * C.H-T- C.N 2 - 
C_iN* — HC.,N— C X H.;- C,N 2 - 

c 2 n*- h 2 — h 2 cn*- c - 

HCO+4- C 2 H 2 — C..HT- CO - 
HCO+4- C 2 H 6 — C_,H7- CO - 
HCO+4- NH:,— NH7- CO - 
HCO+ + hcn — H 2 CX*- CO - 
0+4- No— NO*- N - 
0+4- CH 3 NH 2 — C x H,,Nr- o - 
0+4- CH;,NH 2 — C X H V N7- - 
0+4- CH 3 NH 2 - C X H.,N7- - 

0+4 - ch 4 — ch;- o - 
0+4 - ch 4 — ch;- - 
0+4 - c 2 h 6 — c 2 H 7- - 
0+4 - c 2 h 6 — c 2 Hr- - 
0+4 - h 2 — ohp - H - 

CO++ H — H+ — CO - 

CO++ h 2 — hco*- h + 
co++ ch 4 — ch 7- co - 
CO+4- ch 4 — hco*- ch 3 - 
CO++ ch 4 — c x H7- h -p 
CO+ + c 2 h 2 — CoH*- CO - 
CO+4- NH 3 — NH7- CO -+- 
CO++ O — 0+- CO - 
CO+4- HCN — HCN*- CO - 
CO+4- HC;,N— C x H„N+4- CO - 


1.240 

X 

10“ 

9 

3.800 

X 

10“ 

10 

1.000 

X 

10“ 

• LO 

1.200 

X 

10" 

9 

4.000 

X 

10- 

10 

3.000 

X 

10“ 

10 

1.780 

X 

10“ 

9 

3.150 

X 

10“ 

10 

3.120 

X 

10“ 

9 

7.800 

X 

10“ 

10 

1.190 

X 

10“ 

9 

1.400 

X 

10“ 

10 

7.000 

X 

10“ 

1 1 

6.300 

X 

10“ 

10 

6.500 

X 

10“ 

10 

4.200 

X 

10“ 

10 

7.000 

X 

10' 

11 

2.100 

X 

10" 

10 

1.470 

X 

10“ 

9 

1.280 

X 

10“ 

10 

1.200 

X 

10“ 

10 

3.000 

X 

10“ 

10 

3.600 

X 

10- 

10 

1.200 

X 

10" 

10 

3.000 

X 

10“ 

10 

1.900 

X 

10- 

9 

1.900 

X 

10“ 

10 

4.100 

< 

10- 

9 

3.500 

X 

10" 


8.100 

X 

10- 

O 1 

1.300 

X 

10“ 

'■* 

1.200 

X 

10“ 

0 1 

1.900 

X 

10' 

J 


3.700 ■< 10" y 

1.200 x 10- 1 -(T n ,300)- n "exp(-0.00/T n ) 
1.260 x lO -1 " 

1.660 x lO -9 
3.150 x 10* 10 
8.900 x 10* 10 

1.100 x 10" 10 
1.330 x 10* y 

5.700 x 10“ 10 
1.200 x lO -9 

7.500 x 10 -10 

1.300 x lO -9 
9.790 x 10“ 10 
2.990 x 10“ 10 
8.160 x 10“ 1 1 

4.100 x 10- 10 

1.500 x 10" 9 

1.400 x 10- 10 

3.400 x 10- 9 

2.300 x 10“ 9 
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H+ + E — H + + 

Ht + E — ► H + H - 
H+ + E — H 2 + H - 
Hj+ E — H + H - H 
C++ E — * C + + 

CH++ E — C + H - 
CHt + E — > CH + H — 
CH++ E — CH 2 + H - 
CH++ E -» CH 3 + H - 
CH++ E — CH 2 + H - H 
CH++ E — CH 4 -r H - 
CH++ E — CH 3 + H 2 - 
CH+ + E — CH 3 + H - H 
CH++ E — CH 2 + Ho- H 
Ct+ E — ► C + C — 

C 2 H++ E — ♦ CH + C — 

c 2 h++ e — c 2 h+ h - 

C 2 Ht+ E — CH - CH - 
C 2 H++ E — C 2 H 2 + H - 
C 2 H++ E — CH 2 + CH - 
C 2 H++ E -» CaHs-r H - 
C 2 H++ E — CH 2 - CH 2 -f 
C 2 H++ E — C 2 H 4 - H - 
C 2 H+-t E — CH 3 - CH 2 - 
C 2 H^+ E — CH 4 - ch 2 - 
C X H+ - E — CH 4 - H>- 
N T + E — N + - 
HCN + - E — * CN — H - 
H 2 CN+- E — HCN - H - 
H 2 CN + - E — HNC - H - 
HoCN+- E — CN - H - H 
H 2 CN+- E — CH 2 - N - 
H 2 CN + - E — H 2 — CN - 

C X H V N+ + E — * 

CN+'-u E — C + N - 
C.N+ + E — C + CN - 
C-_»N++ E — C 2 -r N - 
HCO++ E — > H + CO - 


4.000 x 

2.300 x 

4.400 x 
5.600 x 

4.000 x 

1.500 x 

2.500 x 

3.500 x 
1.750 x 
1.750 x 

3.300 x 

1.100 x 
3.300 x 

3.300 x 

3.000 x 

6.500 x 
1.900 x 

8.100 x 

2.300 x 
2.300 x 
1.750 x 
1.750 x 
3.700 x 
3.700 x 
3.700 x 

3.000 x 

1.000 x 

5.000 x 
5.530 x 
5.530 x 
5.530 x 
9.240 x 
9.240 x 

5.000 x 

3.000 x 

1.500 x 
1.500 x 

2.400 x 


C,H + + E — C 4 H 2 - H - 6.200 x 

NH++ E — NH 3 — H - 3.800 x 


0.0 0.0 0.0 


NH++ E — NH 2 + H - 

1.650 x 

NH++ E — ► NH + H 2 - 

1.650 x 

NH++ E — NH + H - 

3.000 x 

NH + + E-N + H* 

4.300 x 

N++ E -» N( 2 £>)+ N - 

1.940 x 

N++ E — N( 2 Z9)+ N( 2 D)+ 

2.600 x 

N 2 H++ E — N 2 + H - 

1.700 x 

CO++ E — C + 0 - 

1.600 x 

C 2 H++ E — C 2 H 6 t h - 

1.000 x 

C 2 H++ E -* C 2 H 5 t h 2 - 

1.000 x 

c 2 h++ E -> c 2 h 4 + h 2 - h 

1.000 X 

H + H — ♦ H 2 + + 

9.030 x 


10- 12 (T e /300)-- r>O0 

10- 7 (T e /300)--" )o 

10- 8 (T,/300)-- 5no 

10 -8 (T p /300) - ' ,o ° 

10 -12 (T f /300) ooo ° 

10- 7 (T p /300)- ,2 ° 

10- 7 (T p /300)- 500 

10- 7 (T,/300)- 500 

10 _7 (T e /300)~ r>o ° 

10- 7 (T e /300)-- 5O ° 

10- 7 (T e /300)-- 27 ° 

10- 7 (T e /300)-' 27 ° 

10- 7 (TJ300)-- 270 

10- 7 (T e /300)-- 27 ° 

10- 7 (T e /300)-- 5o ° 

10- 8 (T e /300)- 500 

10- 7 (T e /300)- r,o ° 

10- 8 (T e /300)- 5O ° 

10 -7 (T e /300)-- 5O ° 

10- 7 (T e /300)- 5O ° 

10- 7 (T e /300)-- 5O ° 

10- 7 (T e /300)- 5O ° 

10- 7 (T e /300)- 5O ° 

10- 7 (T,/300)- >500 

10 -7 (T e /300) - ' r,o ° 

10- 7 (T f ./300)-' 5l, ° 

10-' 2 (T,/300)°- 000 

io- 7 (r„/30o)-- 500 

1 0 - ' ( 77. 300)- 

10“ 8 (T f ;300)- r> °° 

10- 8 (T,/300)- 500 

1 0 _ 8 ( T, / 3 0 0 j ~ • 5 1 w 

1 0 ~ 8 ( T r / 3 0 0 ) “ ' 500 

10 -7 (T,/300j~ - ’ ,o ° 

10- 7 (T P /300)- SO ° 

10- 7 (T,/300)-- 500 

10- 7 (T,/300j-- 5O ° 

io- 7 (r,/30o)-- li9 ° 

1 0 “ 7 ( T,. / 300 j ~ ' 500 
io- 7 (r r /30o)-- 5 °° 


10- 7 (T f /300)-- 500 

io- 7 (r,/3oo)- 500 

io- 7 (r P /3oo)- 500 

10- 8 (T f /300)- 5no 

io- 7 (r f /3oo)- :,,JO 

10- 8 (T f /300)-- 39 ° 

io- 7 (r c /30o)- r> °° 

10 _7 (T f /300) --500 
10- 7 (T e /300)- 5O ° 

10- 7 (T f /300)- 5O ° 

10 -7 (T e /300) _r ’ o ° 

10 _;i3 (T„/300) _ l ' 30 exp(— 0.00/T„) [M] 
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N( 2 L>)4 N 2 — N - No- 
N( 2 Z?)-<- CH 4 — NH + CH 3 - 
N( 2 £>)4 C 2 H 2 — CHCN - H - 
N ( 2 D)+ C 2 H 4 — NH 4 C 2 H 3 t 
N( 2 Z>)+ NH 3 — NH - NHo- 
N( 2 £>)+ H 2 — NH 4 H - 
N( 2 £>)+ CO — N 4 CO - 
N( 2 Z>)4 O — N 4 O 4 
N( 2 £>j+ NH — N 2 + H - 
N + CH 2 — HCN 4 H - 
N 4 CH 3 — H 2 CN+ H - 
N + CH 3 — HCN 4 H - H 
N + Ho— NH 4 H - 
N + H 2 CN— NH - HCN - 
N + CH — CN - H 4 
N 4 Co— CN + C- 
NH 4 H — N 4 H 2 + 

NH + N — N 2 + H - 
NH + NH — No+ Ho- 
CN - CH 4 — HCN 4 CH 3 - 
CN 4- CoH 6 — HCN -r C 2 H 3 - 
CN 4 C 2 H 2 — HC 3 N+ H - 
CN 4- C 2 H 4 — - H - 
CN 4 HCN — C 2 No- H - 

CN 4 CH 3 — CH 3 CN 

CN 4 CHo— CH - HCN - 
CN - CH,— HCN - CHj- 
N - CN — No- C - 
NH - CN — HCN - N - 
Ho 4 CN - HCN - H - 
HCN - H — HoCN— - 

HCN - CoH— HC : ,N- H - 
HoCN- H — HCN - Ho- 
CHCN — N — CoNo- H - 
CHCN - CHCN — C.(No— Ho- 
C;,N4 CH.,— HC : ,N- CH,- 
C 3 N4 CoH,-,— EIC t N-t- CoH-,— 
HC 3 N- H — H 2 C 3 N- - 

HC ; ,N- CN — C,No— H - 
H 2 C 3 N- H — CoHo- HCN - 
C 2 N 2 - H — HCoNo-r - 

HC 2 N 2 - H — HCN 4 HCN - 
CH 4 CH 4 — C 2 H 4 t H - 
CH 4 O — H 4 CO - 
CH 4 O — OH 4 C - 
CH 4 Ho— CH 3 4 4 
'CHj-r No— CH 2 T No- 
*CH 2 4 Ho— CH 2 4 Ho- 
'CH 2+ Ho— CH 3 4 H 4 
'CHo4 CH 4 — CH 2 4 CH 4 ^ 


9.400 x 10" l4 (T n /300)°' ooo exp(— 510.00/T n ) 

1.500 x 10“ 12 (T n /'300) u ' 5Oo exp( — 0.00/Tn) 

1.000 x 10 -11 (T„/300)°‘ 5OO exp(— 0.00/T n ) 

6.000 x 10- n (T n /300)°- 5OO exp(-0.00/T n ) 

5.000 x 10“ 1 '(T„ /300) u,5O °exp(— 0.00/Tn) 

1.700 x 10“ l2 (T n /300)°‘ 5Oo exp(— 0.00/T n ) 

1.700 x 10" 12 

6.900 x 10“ 13 (T n /300) U-5Oo exp(— 0.00/T n ) 

1.000 x 10" 10 

5.000 x 10“ 1 l exp( — 250.00/T„) 

4.300 x 10“ 10 exp(— 420.00/T n ) 

8.600 x 10“ n 

5.200 x 10“ l0 exp(— 16000. 00/T n ) 

3.900 x 10“ u 

2.100 x 10“ ll (T n /300)°' 5OO exp(— 0.00/T n ) 

5.000 x 10“ ll (T n /300)°' 5OO exp(— 0.00/T n ) 

8.200 x 1 0“ 11 (T„/300) a68o exp(— 950.00 /T n ) 

1.900 x 10“ lo (T e /300)°' 5O ° 

1.000 x 10 _33 (T„/300)°' oo exp(— 0.00/T n ) [M] 

7.300 x 10 _13 (T n /300) 2 ' 3OO exp( — 16.00/T n ) 

5.600 x 10“ ^(Tn/SOOj'^^exp^— 500.00/T n ) 

2.900 x 10-‘°(T„/300)- 33O exp(-0.00/T n ) 

2.500 x 10-‘°(T„/300)- ‘ MO exp(-0.00/T n ) 

1.800 x 10“ 11 

1.000 x 10“ 16 

5.300 x 10“ l2 (T„ /300V u)Oo exp( — 2500.00/T,,) 

9.210 x 10“ 12 (T n /300) o ' 70o exp(— 1500.00/T n ) 

1.000 x 10“ lo (T n/ 300)°' O0o exp(— 50.00 /T n ) 

2.940 x 10“ l2 (T n /300) IJ ' 50o exp( — 1000. 00/Tp) 

2.600 x 10“ 1 '(T n /300) o,ooo exp( — 2070.00/Tn) 

k n = 7.110 x 1 0 “ 3l '( 3 0 0 /T rl ) 2 expi — 1 200 . 0 0 / T n ) [M] 
Jfc.* = 9.200 x 10“ 12 exp(- 1200.00 /T n )>ll 

2.200 x 10“ 12 

7.000 x 10“ 11 

1.000 x 10“ 12 

5.000 x 10“ n 

5.000 < 10" " 

1.000 x 10“" 

i fc 0 = 7.110 x 10- ;iO (300/T n ) 2 exp(-1200.00/T„):Ml 
Jt^ = 9.200 x 10“ l2 exp( — 1200.00 /T n )[M] 

1.000 x 10“" 

1.500 x 10“" 

ifco = 7.110 x 10“ 3O (300/T n ) 2 exp^-0.00/T n )[Ml 
k x = 1.500 x 10“ l5 exp(— 0.00/T n )[M] 

1.700 x 10“ 13 exp( — 1 10.00/T n ) 

1.000 x 10“ l ° 

4.750 x 10“ ll (T„/300)°' 5Oo exp( — 0.00/T n ) 

4.750 x 10“ 1 1 (T„ /300) o5OO exp( — 0.00 /T n ) 

2.500 x 10“ 3O (T n /300) ,, oo exp( — 0.00/T n ) [M] 

7.900 x 10- 12 

1.000 x 10“ 12 

7.000 x 10“ 12 

1.600 x 10" 12 
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‘ch 2 + ch.,-> ch 3 + ch 3 - 
ch 2 + ch 2 -> c 2 h 2 + H 2 -i- 

CH 2 + C 2 H 2 — C3H4 + - 
CH 2 + H — CH - Ho- 

ch 3 + ch 3 — c 2 h 6 + - 

CH 3 + H — CH^T + 

C 4 H 2 + H — C 4 H 3 + - 

C 4 H 3 + H — C 4 H 2 + H 2 + 
C 2 + H 2 — C 2 Ht H + 

C 2 + CH.,— C 2 Ht ch 3 + 
c 2 h+ H 2 — C 2 H 2 -f h + 
c 2 h+ ch 4 — c 2 h 2 + ch 3 + 
c 2 h+ c 2 h 2 — c 4 h 2 + h + 

C 2 H 2 + H — ♦ C 2 H 3 t — 

c 2 h 4 + h — C2H5+ - 
c 2 h 3 + h — c 2 h 2 + h 2 - 

C 2 H 5 + H — CH 3 + CH;,-r 
C 4- H 4 — CH 1 "— H 2 - 
C + HJ— CH* — H - 
C + CH* — C.H*- H>- 
C + CHr- CH“~ CH.,- 
C + H 2 CN* — CHT - CN - 
C -i- CN* — C*- CN - 
C 4- CHT — C 2 H“- H - 
C - NHT— * CH*- NH - 
C + HCN 1 " — » CH* - CN - 
C + C)H*— * C r II*- H - 
C 4 - N 2 H* — CH*- N 2 4- 
C -r NHJ — * CHT— NH - 
C + C 2 HT — C r H,*- H -4 
C + CHr — CH*- CH t - 
C + NH — CN - H - 

c - nh 2 — + H - 
c + h 2 — ch 2 

C + N 2 — CNN - - 
C + CO — + - 

C + N — CN 

C + C — ► C 2 t — 

C + H — CH — - 
C + O — CO - - 
C + CH — C 2 + H + 

C + CNN — CN +CN + 
CNN + N — CN + N 2 + 
CNN + H — CH 4. N 2 + 
CNN + H — CN - NH + 
CNN + C 2 H 2 — - + 

CNN + C 2 H 4 — - + 


1.900 x 10“ 12 
5.300 x 10"“ 

k n = 7.600 x 10“ 25 (300/T n ) 2 exp(— 0.00/T n )[M] 
ifc. ^ = 5.900 x 10“ !2 exp( — 0.00/T n )[M] 

4.700 x 10“ lo (T„/300) aooo exp( — 370.00/T n ) 

ko = 1.300 x 10“ 23 (300/T n ) 2 exp(— 0.00/T n )[Ml 
k-c = 5.500 x 10“ 1 *exp(— 0.00/T n )[M] 
kr\ = 1.700 x 10“ 27 (300/T„) 2 exp(— 0.00/T n )[Mi 
= 1.500 x 10“ 10 exp(— 0.00/T n )[M] 
k 0 = 1.390 x 10“ lo (300/T n ) 2 exp(-1.18/T„)iM] 

Jfc^ = 1.000 x 10“ 12 exp(-1.18/T n )[M] 

1.200 x 10““ 

1.400 x 10“ 12 

1.900 x 10“ 11 

1.900 x 10“ n exp(— 1450.00/T n ) 

2.800 x 10“ 12 exp(— 250.00/T n ) 

6.700 x 10“ 11 (T n /300)°' 5OO exp(— 0.00/T„) 

fco = 7.110 x 10“ 3O (300/T„) 2 exp( — 1200.00/T n )[Ml 
k x = 9.200 x 10“ 12 exp( — 1200.00/T n )[M] 

An= 1.220 x 10“ 28 (300/T„) 2 exp(-1040.00/T„)[Ml 
k x = 3.700 x 10“ u exp( — 1040.00 /T n )[Ml 
1.500 x 10- 11 

1.900 x 10“ 10 exp(— 440.00/T n ) 

2.000 x 10“ 9 

2.400 x 10“ 9 

1.200 x 10“ 9 

1.200 x 10“ 9 

1.200 x 10“ 9 

1.100 x 10“ 9 

1.200 x 10“ 9 

1.200 x 10“ 9 

1.100 x 10“ 9 

1.100 x 10“" 

1.100 x 10“ 9 

1.000 x 10““ 

1.100 x 10“ 9 

1.200 x 10“ 9 

2.000 x 10““ 

2.000 x 10“ u (T n /300)°- 5OO exp(-0.00/T n ) 

1.420 x 10“ 31 (T n /300)“ 2OO exp(— 0.00/T n ) [Ml 

6.200 x 10- 33 (T„/300)- 2 - oo exp(-0.00/T n ) [Ml 
1.260 x 10- ;u (T n /300)- 2 - oo exp(-0.00/T n ) [Ml 

1.000 x 10“ 17 (T„/300)°‘ ooo exp(— 0.00/T n ) 

1.000 x 10“ 17 (T n /300) oooo exp( — 0.00/T n ) 

1.000 x 10“ l7 (T„/300) o ooo exp( — 0.00/T n ) 

2.100 x 10“ 19 (T„/300)°' ooo exp( — 0.00/T n ) 

1.730 x 10- u (T„/300)°- 5OO exp(-0.00/T n ) 

1.000 x 10“ 10 (T n /300)°' ooo exp( -0.00 /T n ) 

1.000 x 10“ lo (T„/300) O-OOO exp( — 0.00/T n ) 

5.000 x 10“ ll (T n /300) o ooo exp(-0.00/T n ) 

5.000 x 10“ 11 (T„/300) o ooo exp( -0.00/Tn) 

1.000 x 10“ ,o (T„/300) o ' ooo exp(— 0.00/T n ) 

1.000 x 10“ lo (T n /300) o ooo exp(— 0.00/T n ) 
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H 2 + NH,-> NH'i- H - 
c + nh 2 -> + H - 

O + NH 2 — -r H - 
O + NHo— NH - OH - 


6.320 x 10- u (T„/300)°- 176 exp(— 0.00/T n ) 
2.000 x 10“ 1 1 (T„ /300)°' 5OO exp(— 0.00/T n ) 
3.200 x 10- 1 -(T n /300)°- ooo exp(-0.00/T rl ) 
3.500 x 10“ 12 (T„ /300)°' 5OO exp(— 0.00/T n ) 
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77>/ e _y Selected Proton Affinities 0 


.Neutral 

Ion 

Proton 

Species 

Produced 

Affinity (eV) 

H 

H+ 

2.69 

X 

XH + 

3.37 

h 2 

HJ 

4.39 

CO 

COH+ 

4.74 


XoH + 

5.13 

ch 3 

ch 4 + 

5.62 

ch 4 

ch 5 + 

5.72 

c 2 n 

c 2 nh+ 

5.75 

OH 

H,O t 

6.16 

c 2 h 6 

C 2 Ht 

6.23 

c 2 h 5 

C 2 Hj 

6.42 

c 3 h 8 

C 3 H 0 

6.50 

NH 

NHL 

6.62 

c,h 2 

C 2 H 3 - 

6.69 

c 3 h 7 

c 3 h s ^ 

7.03 

c-,h 4 

C.HT 

7.05 

h 2 o 

H 3 0- 

7.24 

HCX 

HCXH" 

7.46 

c 2 h 3 

C>HT 

-± 

7.56 

hx 3 

H 2 xr 

7.76 

ch 3 cch 

CHiCCH.7 

7.78 

c-C 3 Hc 

i c-C 3 Hg')H~ 

7.80 

c 4 ho 

C 4 H' 

7.79 

c 3 n 2 

C 3 XoH + 

7.81 

c 2 h 

C 2 H7 

7.95 

xh 2 

XHf 

8.13 

HNC 

HCXH + 

8.14 

CH 3 CN 

CH 3 CNK~ 

8.18 

ch 2 chcn 

CHoCHCXH* 

8.23 

ch 2 

CHJ 

8.56 

nh 3 

XHt 

8.85 

x,h 4 

X 2 Hj 

S.87 

ch 3 nnch 3 

ch 3 xhxch 3 + 

8.97 

ch 3 chnh 

ch 3 chnh? 

9.15 

ch 2 chnh 2 

CH 3 CHXHJ 

9.37 

ch 2 nch 3 

CH 9 XHCH 3 

9.41 


a Computed from data taken from Lias et al. 11988] 
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